Cyclotron motion in the vicinity of a Lifshitz transition in graphite 
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Graphite, a model (semi)metal with trigonally warped bands, is investigated with magneto- 
absorption experiment and viewed as an electronic system in the vicinity of the Lifshitz transition. 
A characteristic pattern of up to twenty cyclotron resonance harmonics has been observed. This 
large number of resonances, their relative strengths and characteristic shapes trace the universal 
properties of the electronic states near a separatrix in the momentum space. Quantum-mechanical 
perturbative methods with respect to the trigonal warping term hardly describe the data which are, 
on the other hand, fairly well reproduced within a quasi-classical approach and conventional band 
structure model. Trigonal symmetry is preserved in graphite in contrast to a similar system, bilayer 
graphene. 



Lifshitz transition [lj (also known as electronic topo- 
logical transition) is a change in the Fermi surface topol- 
ogy occurring upon a continuous change of some exter- 
nal parameter, such as pressure Q, magnetic field Q 
or, most naturally, doping This transition does not 
involve a symmetry breaking, alike conventional phase 
transitions of the Landau type, but still leads to observ- 
able singularities in thermodynamics, electron transport, 
sound propagation, and magnetic response of metals Q . 
Saddle points in electronic dispersion, often apparent in 
complex metals, have only recently been visualized with 
the spectroscopy method of angle-resolved photoemis- 
sion [6[. In this Letter, we show how the proximity to a 
Lifshitz transition manifests itself in cyclotron resonance 
(CR) absorption experiments on graphite, a model sys- 
tem with saddle points due to the trigonal warping of 
electronic bands [7|. 

Classically, CR can be understood from the equation 
of motion for an electron in a magnetic field B 



dp/dt = (e/c) [v x B] 



(1) 



where p = ftk is the electron quasi-momentum, e = — |e| 
the electron charge, and v = de(p)/dp is the electron 
velocity, determined by the dispersion e(p). Since both 
the energy e) and the momentum component p z along B 
are conserved, the motion occurs along cyclotron or- 
bits in the {p x ,Py) plane, determined by the condition 
t{Px,Py,Pz) = const. This motion is periodic, and its 
period, 2it/uj c , being proportional to the cyclotron mass, 
defines the cyclotron frequency u> c = uj c (e,p z ). When 
an electric field, oscillating at frequency u, is applied, 
the electron can absorb energy. Absorption becomes res- 
onant when the perturbation frequency w matches the 
cyclotron frequency tu c or its integer multiple. 

In good metals, the incoming radiation is efficiently 
screened and penetrates the sample only within a thin 
skin layer. CR absorption is then a surface effect, ob- 
served mainly when the magnetic field is parallel to the 



surface 0, [T3|. This makes CR for good metals a less 
efficient tool for probing the Fermi surface, as compared 
to other methods, such as, e. g., de Haas-van Alphen ef- 
fect. The resonant absorption is also often smeared by 
the dependence of u> c on p z , which is an additional dis- 
advantage. 

We have applied the CR absorption technique to study 
the cyclotron motion in the vicinity of the Lifshitz transi- 
tion in graphite. The low-temperature in-plane conduc- 
tivity of this material is relatively low, a ~ 10 7 — 10 s (il ■ 
m) , and it quickly decreases upon the application of 
a magnetic field @, EU]. The skin depth thus reaches 
tens of nanometers and greatly exceeds the spacing be- 
tween adjacent graphene layers. Moreover, graphite is a 
highly anisotropic crystal with rather flat electronic dis- 
persion in the z direction (perpendicular to the layers). 
It appears as a suitable material for CR studies of the 
electronic system near the Lifshitz transition driven by 
the trigonal warping of electronic bands. 

CR absorption was measured using the setup rou- 
tinely applied to high- freq uency electron paramagnetic 
resonance experiments 12j. A flake of natural graphite 



(50 (im thick, area 1mm 2 ) was placed in a Fabry- Perot 
cavity mounted inside a superconducting coil. The mag- 
netic field was applied perpendicular to the graphene lay- 
ers. Linearly polarized microwave radiation from a Gunn 
diode tripled to frequency of 283.2 GHz (1.171 meV) 
was delivered to the sample via quasi-optics waveguides. 
The field-modulation technique was applied to enhance 
the detection sensitivity. The modulation amplitude was 
chosen in a way to maximize the signal but to not distort 
the measured lineshapes. 

A representative experimental spectrum (raw data) is 
shown in Fig. [Ha). This trace represents the response of 
the natural graphite specimen measured as a function of 
the magnetic field at fixed microwave frequency. Because 
of the field modulation technique, it corresponds to the 
derivative of the absorbed power with respect to B. The 
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FIG. 1. (color online) Magneto- 
absorption spectra of natural graphite 
measured at a fixed microwave excita- 
tion energy hw = 1.171 meV and de- 
tected with the help of the field modu- 
lation technique at temperature of 5 K. 
Harmonics of fundamental CR frequency 
ujco = eB/m c o (m c0 = 0.057m ) are 
observed down to fields of 20 mT. 
(a) Derivative of the absorption with re- 
spect to the magnetic field B, as a func- 
tion of B. (b) The same plotted as a 
function of uj/oj c o, so that individual har- 
monics at frequencies of |3fc ± 1|cj c o, k = 
0, ±1, ±2 . . ., are clearly seen, as marked 
by vertical arrows, (c) Absorption as a 
function of uj/ u) c o obtained by the numer- 
ical integration of the curve presented in 
part (a) with respect to B. 



magneto-absorption response of graphite is expected to 
be mostly sensitive to singularities in the electronic joint 
density of states, located at the K and H points of the 
graphite Brillouin zone. A number of the observed reso- 
nances can be easily identified as due to electronic states 
at the K point, along the results of previous similar stud- 



ies [13l - ll8| . Holes at the H point as well as decoupled 
sheets of graphene on the surface of graphite give rise to 
resonances at different spectral range (much lower mag- 
netic fields) mm. 

To a very first approximation, the K point electrons 
of graphite have parabolic dispersion. Their effective 
mass, most frequently reported to be in the range from 
m = 0.057m to 0.060m Q (mo free electron mass), 
fixes the cyclotron frequency at hiu c a k 2 x B[T] meV. 
Then, the broad but still visible resonance at \B\ ~ 0.6 T 
is attributed to the fundamental CR absorption. All 
other observed resonances are higher harmonics of the 
fundamental one. This is evidenced in Fig. [ljb) where 
the spectrum from Fig. QJa) is re-plotted against w/w c0 
(i.e. versus instead of B). frui c o is eventually set at 
2.05xi3[T] meV. In agreement with previous reports [Uj], 
the observed harmonics follow two series: cj w |3fc±l|w c rj, 
where k = 0,±1,±2.... 

The superior quality of the present data (due to 
higher frequencies applied and perhaps better quality of 
graphite specimens) allows us to uncover more and in- 
triguing spectral features. Our key observations, that 
we interpret in the following, are: i) the appearance of 
a large number (up to 20) of CR harmonics ii) an en- 
hanced strength of 3fc + 1 harmonics as compared to the 
strength of the 3k — 1 series at B > (and vice versa at 
B < 0) and finally, iii) a very characteristic, asymmetric 
broadening of the observed resonances, enhanced on the 
low frequency (high-field) sides of the absorption peaks. 
These features are clearly seen in the raw data and also 
in Fig. QJc) in which we reproduce the actual absorption 



spectrum, as derived from the numerical integration, over 
the magnetic field, of the measured (differential) signal. 

The appearance of n-harmonics with n = 3fc ± 1 is 
usually understood as due to breaking of the isotropy of 
the electronic spectrum in the layer plane by the trigonal 
warping. For isotropic bands, only the k = fundamen- 
tal transition is allowed, whereas the n — 3fc±l harmonic 
appears in the |fc|th order of the perturbation theory with 
respect to the trigonal warping term 3, 21 1 . The spec- 
trum in Fig. Q] contains many harmonics which start to 
fall off only at large indices n > 7. Clearly, the pertur- 
bation theory is not applicable to interpret these data. 
Instead, we will use the quasi-classical approximation. 

Eq. ((TJ) can be cast in the Hamiltonian form in 
the phase space (p x ,p y ): dp x /dt = —dH(p x ,Py)/dp y , 
dpy/dt = dl-L(p x ,Py)/dp x , with the Hamiltonian 
T-L{PxiPy) = —(eB/c) e(jp x ,p y ) (we omitted p z , which en- 
ters as a parameter). Generally, classical Hamiltonian 
systems exhibit a very rich behavior. However, they 
share some universal features when the energy e is close 
to that of a saddle point of the Hamiltonian, e = e S p, 
as is well-known in the classical nonlinear physics [22j. 
(i) The cyclotron motion in the vicinity of a saddle point 
is slow and its period diverges logarithmically, u) c {e) — > 
for e — > e sp . (ii) The Fourier spectrum of this motion con- 
tains many harmonics and their number diverges when 
e — > e sp . This second fact provides an obvious hint for 
the interpretation of the experimental data. 

The experimentally probed electronic states are those 
around the Fermi level ef- Thus, the effects discussed 
above are important if ep m e sp . This is the case of 
graphite, as illustrated in Fig. [2] using standard calcu- 
lations based on the Slonczewski-Weiss-McClure (SWM) 
model [23[ in the two-band approximation (see Supple- 
mentary Information). Here we used the standard val- 
ues of the SWM parameters : 7o = 3150 meV, 71 = 
375 meV, 72 = -20 meV, 73 = 315 meV, 74 = 44 meV, 
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FIG. 2. (color online) (a,b) Electronic structure near the 
K point of graphite (k z = 0). Two separatrix lines pass 
through six saddle points. The Fermi level is located about 
6 meV above the upper separatrix. (c) Constant energy 
contours in the (k x ,k y ) plane for k z — for e — —17, 
e = —25 (Fermi level), e = —30.6 (upper separatrix), and 
e — —33 meV. (d) Classical cyclotron frequency hu) c {c, k z = 0) 
at B = 100 mT in the relevant energy interval. uj c van- 
ishes at the saddle point. Open circles show the LL spacing, 
Aej = ei+i — e;, as a function of ej, derived from the SWM 
model. Roughly 1 meV away from the saddle point, the circles 
fall on the classical curve, Ae; w huj c (ei). 



75 = 38 meV, A = —8 meV. The band dispersion has six 
saddle points at two different energies e - sp and eh_ sp , 
which define two separatrices - isoenergetic lines separat- 
ing regions with different topology. Fermi level crossing 
these saddle points would imply the change in the topol- 
ogy of the Fermi surface, which actually corresponds to 
the Lifshitz transition of the neck-collapsing type. The 
Fermi level is close to the upper separatrix, on which we 
focus our attention hereafter, e sp = e c -s P - The single 
electron pocket around the K point at ep > e sp , splits 
into four disconnected pockets when ep goes below e sp . 
Fig. |2fd) shows the classical cyclotron frequency for the 
SWM dispersion at k z =0, which vanishes at e = e sp . 

In the language of quantum mechanics, the k z = en- 
ergy spectrum consists of discrete Landau levels (LLs) 
e;. In the quasiclassical approximation, ei can be found 
from the Bohr-Sommerfeld quantization rule. The n th 
CR harmonic corresponds to the transition over n lev- 
els, nhui c rj e;_|_ n — e;, to the leading order in K. The 



decrease of u) c (e) at e — > e sp corresponds to an accumu- 
lation of LLs. Nevertheless, Ae does not approach zero, 
since the condition of the validity of the quasiclassical 
quantization, \u) c (e + hu c ) — w c (e)| <C w c (e), holds only if 
e is not too close to e sp . LLs always remain discrete, see 
Fig. [3d). 



As we will show later, ep — e sp is about 6 meV, i.e., 
it is five times larger than the microwave frequency, 
Hd = 1.171 meV. Our classical approximation is justified 
in this case. As a matter of fact, the quasi-classical ap- 
proximation works the better, the smaller uj is. However, 
if microwave frequency is too small, the harmonic struc- 
ture will be smeared by broadening of electronic states. 
The optimal frequency, used in the experiment, is thus 
determined by an appropriate compromise between these 
two competing conditions. 

Assuming that the absorbed power is proportional to 
the real part of the conductivity, Re a xx {uj), and calcu- 
lating the latter from the standard kinetic equation [l[ in 
the simplest relaxation time approximation for the colli- 
sion integral (see Supplementary Information), we obtain 



Rea xx {uj) 



e 



,— — oo 



df 

de 



Fm c \v x J 2 dk z de 



h 2 (uj~nuj c ) 2 + r 2 ' 
(2) 

where the k z integration is from — 7r/(2a z ) to n/(2a z ). 
Both the cyclotron frequency, w c , and the cyclotron mass, 
m c = —eB/(cuj c ), depend on e and k z . The basic fre- 
quency w c o, introduced earlier, is w c o = ui c (e = ep,k z = 
0). v„ = v„(e, k z ) is the Fourier harmonic of the elec- 
tron velocity, corresponding to the term oc e~ mUct , de- 
termined from the solution of the unperturbed equation 
of motion, Eq. (fl}. Finally, T accounts for relaxation, 
and /(e) is the Fermi function. 

Even without solving Eq. (JTJ , it is easy to see that the 
triangular symmetry of e(p) in the (p Xl p y ) plane fixes 
v„ = for n = 3k, k = ±1,±2,.... In Fig. QJb), 
the resonances at n — 3k are absent, which demon- 
strates that the triangular symmetry is not broken in 
graphite. This is in contrast with recent reports for 
a bilaycr graphene (24-27|, even though it is formally 
described by the same single-particle Hamiltonian (for 
a fixed k z ) [28| . The Sclulg symmetry fixes v n ^ x to be 
real and v„, y = ±iv n ^ x for n = 3k ± 1, so the peaks at 
n — 3k + 1 and n = 3k — 1 are seen in the opposite circu- 
lar polarizations of the microwave field. This helps us to 
understand the observed difference in the intensities of 
the n = 3k + 1 and n = 3k — 1 series. Indeed, for B > 0, 
when the electron moves along the Fermi surface, shown 
in Fig. Efb), in the overall counterclockwise direction, it 
should be more strongly coupled to the counterclockwise 
polarized radiation. As both circular polarizations are 
equally present in the incoming linearly polarized radi- 
ation, the spectrum is fairly symmetric with respect to 
B -> —B. 

The present studies are restricted to bulk graphite, a 
system with fixed Fermi level but in an apparent prox- 
imity to the Lifshitz transition. An obvious experimen- 
tal challenge would be to trace the CR response when 
changing the Fermi level with respect to the separa- 
trix energy, with an attempt to tune the proximity to 
Lifshitz transition in graphitic structures. This can be 
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FIG. 3. Derivative of the absorption with respect to B, 
as calculated using Eq. ((2j) for different electronic broaden- 
ings r. (a) a constant broadening T — 20 /ieV and (b) V = 

O.lV^BpTmeV. 



in principle envisaged for electrostatically gated bilayer 
graphene [2(| l27ll a nd /or for bulk graphite under hydro- 
static pressure [29|. Importantly, such experiments re- 
quire no degradation of the quality of the sample, which 
likely excludes the experiments on, for example, chemi- 
cally doped structures. 

Besides the large number of harmonics, typical of a 
classical motion near a saddle point, the proximity to 
the Lifshitz transition also leads to some lowering of the 
cyclotron frequency. Indeed, the fundamental cyclotron 
frequency determined from the period in B^ 1 of the spec- 
trum, hu; c o/B = 2.05 meV/T, is slightly lower than its 
parabolic-band limit at k z = 0, 2.24 meV/T. The latter 
value, however, relies on the specific values of the pa- 
rameters of the SWM model. More apparent effects are 
deduced from the analysis of the peak shapes (which are 
determined by the integration over k z ). 

As seen from Fig. [Qc) , each peak has an abrupt cut-off 
on the high-frequency side and a tail on the low-frequency 
side. This contradicts the first intuition, based on the 
well-known fact that the parabolic part of the bands be- 
comes steeper as k z increases from the K point towards 
the H point. This pushes the LLs upwards as k z increases 
and would result in a tail on the high-frequency side 
of each peak in the absorption spectrum A(oj,B) [3(ij |. 
However, the bottom of the conduction band [defined 
as e(p = 0, k z )} and the saddle point e sp (k z ) shift up- 
wards upon increasing fc z , as e(p = 0, k z ) = 2j2 cos k z a z , 
72 < 0. Thus, the Fermi level approaches the saddle 
point as k z moves away from k z = point and uj c (ep,k z ) 
decreases simultaneously, see Fig. [2jd). This provides a 
tail on the low-frequency side of the peaks. Thus, the 
suppression of lo c near the Lifshitz transition is crucial to 
interpret the peak asymmetry. 

The spectrum derived from Eq. ([2]) is shown in Fig. 02 
with ep and T as the only adjustable parameters - the 
parameters of the SWM model were fixed 0- The 



best agreement is obtained for ep = —25 meV. If 
ep = —24 meV, the peaks have no asymmetry, since 
u c (cf = — 24meV, k z ) has a significant upturn on in- 
creasing k z . When ep — — 26meV the fall-off of largc- 
n harmonics is noticeably slower than the experimental 
one. In other words, the closer ep is to e sp , the more 
harmonics are seen in the spectrum. The frequency 
w c (e F = -25meV,A; 2 = 0)/B = 2.03mcV/T agrees 
with the experimental value, 2.05 meV/T. The value 
ep = —25 meV is also in good agreement with the one de- 
termined independently from the charge neutrality con- 
dition (see Supplementary Information), ep = — 24meV. 
A constant value of V = 20/ieV was assumed for the 
curve in Fig. EJa). Apparently, this does not describe 
well the amplitudes of peaks at n = 1,2: the theoreti- 
cal peaks are narrower and thus more intense than the 
experimental ones. Better agreement is obtained under 
assumption that T cx y/B (see Ref. HI and Supplementary 
Information). Notably, the curve in Fig. EJb) with T = 
0.1\/ B[T] meV corresponds to the zero-field relaxation 
rate H/tb=o = 40 /icV. The extracted value of tb=o pro- 
vides the zero- field dc conductivity a = 1.2x 10 8 (i7-m) _1 
(see Supplementary Information). This is fully consistent 
with typical literature data @, EH and implies a mean 
electron free path of 6 /im, which is, notably, comparable 
or even longer than the corresponding values reported for 
strictly 2D graphene-based structures (2(| (3^, HI] . 

To conclude, we have introduced CR experiments as a 
new tool to study Lifshitz transitions. We have shown 
how the proximity to the Lifshitz transition manifests it- 
self in the CR spectrum of a model system, bulk graphite. 
Namely, we have observed a multi-mode response, where 
the basic CR mode is accompanied by many harmonics. 
Using the standard SWM model for the electronic band 
structure of graphite to analyze the data, we have de- 
termined the Fermi energy and estimated the electronic 
broadening. The similarity between the band structure of 
graphite near the K point and that of a bilayer graphene 
logically suggests to probe Lifshitz transition in the latter 
system by CR methods, and to shed more light on the 
currently debated issue of spontaneous symmetry break- 
ing in bilayer graphene [24H27| . 
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SLONCZEWSKI-WEISS-MCCLURE MODEL AND 
ITS TWO-BAND PROJECTION 

The Hamiltonian of the AB-stacked graphite can be 
written as 
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(3) 



where each element is a matrix in the Hilbert space of 
a single layer. H± and H 2 contain the in-plane nearest- 
neighbor coupling with the matrix element 70 as well as 
the diagonal shift A. V\ 2 represents the coupling between 
neighboring layers with matrix elements 71 for the neigh- 
boring atoms, and 73, 74 for the second nearest neighbors. 
V\\ and V 22 correspond to second- nearest-layer coupling 
with matrix elements 72/2 and 75/2. Due to the transla- 
tional invariance in the z direction, the problem can be 
reduced to that of an effective bilayer: 



Hi + Vne 2lk * a - + V n e~ 2ik > c 
2V12 cos k z a z 



2V12 cos k z 



H 2 + V2 2 e 2lk * a * + V". 



t —2ik z a,z 



(4) 



where a z is the distance between the neighboring layers, lar to the layers, —ir/(2a z ) < k z < ir/(2a z ) (note that 
and k z is the wave vector in the direction perpendicu- 
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the period of the structure is 2a z , that is, two layers). 
Expansion of the single-layer Hamiltonian in p x ,p v , the 
in-plane quasi-momentum components counted from the 
H — K — H line, gives 



HkAp) 



where 



«3,4 



3 70a 
2 h ' 
273,4 
7o 



r 2 


vp~ 


— a 4 vp- 


a 3 vp + 


vp+ 


r 5 




—a^vp- 


-a 4 vp^ 


- ri 


r s 


Vp- 


a 3 vp- 


—a 4 vp + 


vp + 




Ti = 


271C, r 2 


= 272C 2 , 





(5) 



c, r 5 = 2 75 c 2 



C = cos k z a z 



a = 1.42A is the distance between the neighboring carbon 
atoms in the same layer, and p± = —iK(d x ± id y ). It 
is convenient to rotate the basis in the space of the 4- 
columns (ipi,ip2,ip3,ipi) T as 



</4 



j>2 + ^3 

V2 



l/>' 2 = ll>l, "03 = "04, ^4 



j'2 - 03 

V2 



so the transformed Hamiltonian becomes 



T 5 + Ti v 4 p + v^p- 

v 4 p + -w 4 p_ r 5 - r x 



(6) 



(7) 



where u 4 = v (1 + a A )/V2, v 4 = v(l - a 4 )/V2. The four 
eigenvalues for each p x ,p y , k z can be found by calculating 
the determinant 

deti^(p) = (l-a 2 ) 2 M 4 - 

- 2a 3 [(1 + aDTx + 2a 4 r 5 ] {vpf cos 3^ p + 

+ [a 2 (T 2 - T 2 ) - 2(1 + a 2 )r 2 r 5 - 4a 4 r 2 r 1 ] {vpf- 

-r 2 (r 2 -r 2 ), 



r. 



r. 



(8) 

e to obtain 



and replacing 
det(H' kz -e). 

To obtain the effective 2x2 Hamiltonian, acting in 
the subspace of the two low-energy bands, we eliminate 
ip'i,ip' 4 from the Schrodinger equation H' k (p)ip' = eip'. 
This gives 



< X2 (P) = 



where 



T 2 a 3 vp^ 
a 3 vp- - v 2 p'i/Ti 



- v 2 p 2 _/r 1 
r 2 



"Keh 



v*p-p+ 




(eh = "4 





Vp+p- 



2Tt 



(9) 



determines the electron-hole asymmetry of the spectrum. 
At energies we arc interested in (a few tens of meV) the 
two-band approximation works well for all k z except the 
immediate vicinity of the H point, tt/2 — \k z \a z < 0.05. 
In the two-band approximation, there is an analytical ex- 
pression for the energies, which is most conveniently writ- 
ten in the polar coordinates, p x = pcosip, p y = psini^: 



(p cos ip, p sin <p) 



± 



>(vp) 4 

r? 



2a 3 %^cos3<^ + aiM 2 . 
1 1 



(10) 



From this expression one can deduce the energies of the 
saddle point in the conduction band, the energies of the 
three "leg" conical points, where the conduction and the 
valence bands touch each other (in addition to the conical 
point at p = with e = T 2 ), and the energy of the saddle 
point in the valence band: 



4(l-2Cefc)' 

eicg = T 2 + 2Q eh a\Ti, 



^h — sp 



air x 



4(1 +2U)' 



(11a) 
(lib) 
(11c) 



Given the dispersion (|10p . one can determine the con- 
centrations of electrons and holes at zero temperature as 
functions of the Fermi energy, 



n e (e F ) 



e+(p)<e F 



nh(ep) = 4 



E (p)>«F 



d 3 P 

(2irhf ' 
d 3 p 

{2-Khf 



(12a) 
(12b) 



where the factor of 4 takes care of valley and spin degen- 
eracy. If the sample is undoped, e F can be determined 
from the neutrality condition n e (e F ) = nh(€F)- Elec- 
tron and hole concentrations obtained from Eqs. (|12ap . 
(|12b[) , are plotted in Fig. |4l together with the simplified 
expression, 



2 7l \e F \d(±e F )) 
n e ,h{£F) = —5 — To ^- z ( ±4 74/7o) 



Z(x) 



Tr 2 a z (3j a/2) 2 

2 / 7T 1 

x [ 4 ~ Vl - ~ 2 



arctan 



(13) 



obtained by neglecting 7 2 , 73 , 75 . 



CONDUCTIVITY FROM THE KINETIC 
EQUATION 

We determine the conductivity from the electronic dis- 
persion following the standard procedure[TJ. First, let us 
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FIG. 4. Electron and hole concentrations as functions of the 
Fermi energy, as given by Eqs. (|12a[) . (| 12b[) (thick solid lines), 
together with the simplified expression (|13[) (dashed lines). 



consider the problem at some fixed value of k z , and find 
the in-plane electric current produced by an external in- 
plane electric field 



E(i) = E w e~ lult + E* e 



(14) 



The in-plane group velocity of the electrons being v = 
de(p)/dp (p is in the xy plane, and we omit the k z ar- 
gument for the moment), the semiclassical equations of 
motion are 



dp x _ eB de(p x ,p y ) 
dt c dp y 

dp y _eB de{p x ,p y ) 
dt ~ c 



dp x 



eE x (t), (15a) 
+ eE y {t). (15b) 



Here the electron charge e < 0. These equations can 
be written in the Hamiltonian form dp x /dt = —dH/dpy, 
dpy/dt = dl-L/dp x , with the Hamiltonian function 

eB 

H{p x ,Py) = — — e(jp x ,py) + e(E v p x - E x p y ). (16) 

In the absence of the perturbing electric field the trajecto- 
ries in the (p x ,p y )-space (cyclotron orbits) coincide with 
the constant energy contours, e(p) = const. The solution 
of the unperturbed problem (E = 0) can be written as 



Thus defined m c is nothing but the cyclotron mass. The 
harmonics of the velocity can be related to those of mo- 
mentum directly from the equations of motion, which 
gives 



v„ x (e) 



m c (e) 



m c (e) 



Pnx(e). (19) 



The symmetry of the cyclotron orbit with respect to ro- 
tations by 2ir/3 fixes 

P3k,x = P3k,y = 0, P3k±l,x = TiP3k±l,y = ^3fc±l, (20) 

for integer fc, where -Pn(e) are real. When Ch- S p < e < 
e _ sp , in addition to the orbit which is C^-symmetric 
around the origin, for each e there are three other orbits 
encircling the three "leg" conical points, each such orbit 
having only one mirror reflection symmetry. These give 
rise to another series solutions with 



Pn,x — P n , Pn,y — ^Qni 



(21) 



where f^,(e), Q' n ( e ) are a l so reai - However, we will be 
mostly concerned with the first solution, Eq. (|20l) . 

In order to separate the oscillating motion, it is con- 
venient to perform the canonical change to action-angle 
variables of the unperturbed Hamiltonian, (p x ,p y ) — > 
(S, (f). The action variable is defined by 



5(e) = - 



eB 



dp x dp y 
2tt 



(22) 



e(p)<e 



and the angle <f> = u c t for the unperturbed motion. The 
inverse transformation, (S,4>) — > (p X iP y ), is given by 
Eq. (jTT)) with the replacement tu c t —> (f> and e — s- e(S), 
the latter determined from Eq. (f2"2")l . As for any canonical 
transformation, the phase volume is preserved: dp x dp y = 
dS d<f>. In the new variables, the Hamiltonian (|16p as- 
sumes the form 



-irufi 



U(S, 0) = / u c (S') dS' + y eE • v n (S) e 

(23) 

The distribution function J r (5', <f>) is defined as the av- 
erage number of particles in the quantum of the phase 
volume, (27rft) 2 , around some given values of S, <j>. It de- 
termines the total current carried by the electrons: 



p(t)= £ PnWe-^W*, p„(e)=pi„(e), (17) j(t) = f ™** H s, , t ) ]T ev„(5)e 



(24) 



where the frequency is determined by the derivative of 
the area, enclosed by the cyclotron orbit, with respect to 
energy: 



w c(e) = tt-, m c {e) = — — 

m c {e)c 2n oe 



dp x dpy. (18) 



e(p)<e 



The distribution function J- satisfies the kinetic equation 
dF dU dT OH dT 



dt dS 



dS 



= St [71, 



(25) 



where the left-hand side is the canonical Liouville op- 
erator, and the right-hand side represents the collision 
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integral. We seek the distribution function in the form 
F{S,<j),t) = T (S) + 6F(S,<f>,t), where T (S) is the 
equilibrium distribution function, corresponding to the 
Fermi-Dirac distribution /(e): 

MS) = MS))= 1 + elt < L s ) -. r] /T > ( 26 ) 

and the correction ST is sought to linear order in E. For 
the collision term we adopt the simplest approximation 
with relaxation time r, so the linearized kinetic equation 
for ST has the form: 
dSF dSF dF \ „ - , 5F 



t)t 



(27) 

For the monochromatic field, Eq. (|14[) . one readily finds 
the oscillating correction ST(t) = SJ r ul e~ lu ' t + ST^e 1 ^, 
with SJ-^, given by 



dJ^o em c v n ■ E w 



OS 



n— — oo 



-iri(fi 



jr — i(cu + ncj c ) 



(28) 



Substituting this expression in Eq. (|24[) and passing from 
integration over S to integration over e, we obtain the 
current, j(t) = j^e - ^* + j*e iwt , with j w given by 



2irh 2 



"i c v„(v* • E^ 



/r — i(w — jiu c ) 



(29) 



Up to now we considered the problem at fixed k z (recall 
that m c ,tu c , v„ depend on k z via cos k z a z ). The obtained 
current should be summed over k z with the help of 



7r/(2a a ) 



E 



2- 



(30) 



-Tr/(2a z ) 



where N z is the number of graphene layers, and thus 
N z a z is the thickness of the sample. Also, the result 
should be multiplied by the spin multiplicity 2, and by 
the valley multiplicity Af = 2. This gives the in-plane 
conductivity per graphene layer, which should be further 
divided by a z to obtain the bulk conductivity, 



(Tij (u) 



2Afe 2 
2Tth 2 a z 



E 



de 



tt/2 

df\ f d{k z a z ) 



de J J 2tt 

tt/2 



(31) 



m c Vi,n(Vj,n)* 

r — i(oj — nto c ) ' 



where T = 1/r. The B limit is recovered by setting 
cj c — > and 



E 



{VjV.: 



J ViVjS(e kz (p) - e) dpx dp y 
I <K e fc* (p) - e)dp x dp y 



(32) 

We evaluate this integral numerically, and plot the corre- 
sponding dc conductivity in zero magnetic field in Fig. [5] 
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FIG. 5. The in-plane dc conductivity of graphite at zero tem- 
perature and zero magnetic field for h/r = 40 fieV, as a func- 
tion of the Fermi energy. 



SELF-CONSISTENT BORN APPROXIMATION 
IN HIGH MAGNETIC FIELDS FOR ARBITRARY 
ELECTRONIC DISPERSION 

Here we show that the self-consistent Born approxi- 
mation, studied in detail for two-dimensional electron gas 
with parabolic dispersion by Ando and Uemura0,0|, can 
be straightforwardly generalized to arbitrary electronic 
dispersion for short-range disorder and strong magnetic 
field. Moreover, the expressions obtained by Ando and 
Uemura for the broadening of the Landau levels and the 
cyclotron resonance lines remain unchanged. 

Consider two-dimensional electrons with a dispersion 
law e(p). The uniform magnetic field B is perpendicular 
to the plane, and is described by the vector potential A. 
In the presence of a static disordered potential V(r), the 
electronic Hamiltonian takes the form 

H = e(— ifiV — eA/ c) + V(r) =H + H dis . (33) 

The random potential is assumed to be Gaussian with 
zero average, V(r) = 0, and is determined by its pair 
correlator, 



V(r) V(r') = W{r - r') = W S(r - r'), 



(34) 



where the overbar denotes the average over the disorder 
realizations. If the disorder is due to random impurities 
with the two-dimensional density n; mp , and each impu- 
rity has a short-range potential iti m p(r), then 

2 



Wo 



"imp 



*rmp 



(r) d 2 r 



(35) 



Let us define the disorder-free Green's function: 
G (r,r';e) = (r|(e- J ff )- 1 |r') = 

_ ^ V/a(r)V>r«(r') _ ^ Vi(r,r>) (36) 
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Here I labels the Landau levels with energies e;, and a 
labels states on the same Landau level with wave func- 
tions ipia{ r ) (eigenfunctions of Hq). For example, in the 
Landau gauge, A x = —By, they can be labeled by the 
momentum p x : 



^Ivxfav) = eW:c/h Xi{y + c Px/eB), 



(37a) 



where the functions Xl{y) an d Landau level energies e/ 
are found from the following equation: 



e(eBy/c,-ihd y )xi(y) = Qx(y)- 



(37b) 



The choice of basis in the degenerate manifold of each 
Landau level, as well as the specific form of the wave func- 
tions will not be important. We will only need the follow- 
ing general properties: (i) the projector kernel Vi(r,r') 
depends only on the relative coordinate r — r' (and so 
does the Green's function Go), and (ii) the degeneracy of 
each Landau level is L x L y / (2tt£ 2 3 ), where L x and L y are 
the dimensions of the sample in the x and y directions, 
respectively, and £ 2 B = hc/(eB). Thus, 



1 



2<- 



(38) 



In the presence of disorder, one can define the full 
Green's function for each disorder realization, 



G(r,r';e) = (r|(e- J ff)- 1 |r')=X] 



^(r)^(r') 



(39) 



(</> s and E s being the wave functions and energies of the 
exact eigenstates of H), as well as its average over the 
realization, G(r, r'; e), from which several physical quan- 
tities may be obtained, as will be seen below. In the so- 
called self-consistent Born approximation (SCBA), which 
can be represented diagrammatically using the standard 
rulesQ as shown in Fig. [BJb), one obtains the following 
system of closed equations for the average Green's func- 
tion G and the self-energy S jij: 



G(r,r') = G (r,r' 



J G (r, ri)E(n, r 2 )G(r 2 , r') dr 1( ir 2 



E(r,r') = W(r-r')G(r,r'). 



(40a) 
(40b) 



The e argument, common for all functions, has been omit- 
ted for the sake of compactness. SCBA in strong mag- 
netic fields has been analyzed in Ref. [U for the parabolic 
dispersion, e(p) = p 2 /(2m), and was shown to be valid 
for high Landau levels [5j, in agreement with the general 
rule[4|: SCBA is valid when the electron motion between 
successive scattering events is quasiclassical (i. e., con- 
tains many de Broglie wavelengths). Even when this 
condition is fulfilled, the SCBA misses effects related to 
coherent multiple-impurity scattering, but such effects 
are beyond the scope of the present work. 



(a) 



r' 



r' 



r' 



G (r,r';e) 
G(r,r';e) 
W(r,r') 



(b) 



+ 



(c) 



e+co 



+ 



FIG. 6. (a) Notations for different elements, constituting 
the diagrams, (b) diagrammatic representation of the SCBA 
equations (I40a[) . (|40h>[) . (c) SCBA diagrams for the disorder- 
averaged conductivity. SCBA, sometimes also called the non- 
crossing approximation, corresponds to neglecting all dia- 
grams with crossing dashed lines. 



Below we show that for the <5-correlated disorder, 
Eq. (|3"4")l , and for strong magnetic fields, the SCBA equa- 
tions (|40a|) . (|40b[) are straightforwardly solved for any 
dispersion law e(p). Indeed, since W(r — r') = WoS(r — 
r'), and since G(r,r; e) does not depend on r, Eq. (|40bj) 
gives E(r, r'; e) = £(e) <5(r - r'). Thus, E(r, r'; e) has no 
matrix elements between different Landau levels, so 



G(r,r';e) = ^ 



e-Q-E(e)' 



(41) 



For E(e) we have the self-consistency equation, which can 
be written as 



£(e) 



E 7 



7l 



e*-E(e) : 



7l 



W 
2^1 



(42) 



where we used Eq. (|38|) . It is convenient to express 
the factor 7^ in terms of the cyclotron frequency uj c = 
(e; + i — ei)/h and the zero-field scattering rate 1/tb=o, 
assuming that both vary weakly with energy on the scale 
of hw c . The self-consistent Born approximation at zero 
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field gives 



Tb=o 



= 2ttNWqv, 



(43) 



where TV is the number of valleys (provided that the 
disorder is sufficiently short-range to induce efficient in- 
tervalley scattering), and v is the density of states per 
unit area per spin projection and per valley at zero field. 
When the energy dependence of v is weak, one can write 



(44) 



Hu c 2tt£ 2 b ' 



since each Landau level has l/(2ir£ B ) states per unit area, 
and Hw c is the separation between the Landau levels. 
This gives 



1% = 



t b=0 2tt7V" 



(45) 



In the limit of strong fields, lo c tb=q 3> 1, we have 
1b "C £; — £;'■ Let us focus on some Landau level I and 
on energies E close to e/. First, let us consider terms in 
the sum in Eq. (|42p corresponding to levels different from 
the chosen level I. They result in a small overall shift of 
the level I: 



ei + £(ej) w e t + 



E 



Q - ei> 



(46) 



In the following we will neglect the difference e; — e; ~ 
j B /(luo c ), which is much smaller than the Landau level 
broadening, as we will sec shortly. The term in the sum 
in Eq. (|42|) . corresponding to the same level I, should 
be treated exactly; it results in a quadratic equation 
for S(e). As a result, we can write 



E(e) 



e-ei 



IB 



E 

i'^i 



e - ev 



(47) 



for \e - ei\ < hw c . 

The knowledge of S(e), Eq. gTJ), and of G(r,r';e), 
Eq. (|4T|) . enables us to find the disorder- averaged density 
of states per unit area per spin projection and per valley: 



1 



v{e) = -ImG(r,r;e-iO+) 

7T 



ttjb 



c 

2n E 



(48b) 



Thus, each Landau level, which in the absence of disorder 
is discrete and inifinitely degenerate, in the presence of 
disorder is broadened into a semicircle. 

To describe the cyclotron resonance (that is, inter- 
Landau- level absorption), we start from the expression 
for two-dimensional conductivity in terms of the exact 



wave functions <f> s (r) and energies E s , following directly 
from the Kubo formula: 



( 2) 27V ie 2 



E/s 



lJ y "' ujL x L y 

2Mie 2 x . (f.-UMWi^Hjls) 



y 

u)L x L y *—f huj — [E s 



2Mihe 2 



E 



U.-U)WiW)W\*j\*) 



L X L V (E s , - E s )(hu - E s , + E s + iO+Y 

(49) 

Here we used the notations 



de(p) 



dp 



8 2 e(p) 



p— y— ihW 



dpidpj 



(50) 



-ihV 



and the last equality in Eq. (|49| was obtained using the 
commutation relations 



H ,Xj 



(51) 



Introducing the notation 

G R - A (r,r';e) = G(r,r';e + iO + )-G(r,r';e-iO+), (52) 

we can write the dissipative part of the conductivity, av- 
eraged over the disorder, as 



Rea x 2 x \co) 



27Ve 2 

47TW 



de [f(e + fin) - /(e)] 



d 2 vd 2 v' 

L X Ly 



x v x G R - A (r, r'; e + fiu) v' x G R - A (r', r; e). 

(53) 

This expression involves the average of products of two 
Green's functions, which in SCBA reduces to summation 
of the ladder series, shown in Fig. IUc). Consider the 
second term in this series, containing one impurity line. 
It contains the spatial integral 

where ri is the position of the impurity and 

Jete' ± ( r i)= / d 2 r G(n, r; e ± i0+) v x G(r, n; e ± i0+). 

(54) 

Since the average Green's functions depend only on the 
relative coordinate, J^i^ (ri ) does not depend on ri, and 
thus it can be written as 



Jff= I ^ 7 I1 G(r u r;e±iO+)v x G(r,v 1 ;e±iO+). 

J J-ixJ-iy 

(55) 

Recalling Eq. (|4"Tj) for the Green's function, we note that 
because of the integration over ri, the contribution to 
transitions between different Landau levels I and V in- 
volves the product of the projectors V{Pv = 0. Thus, 
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Explicitly, from Eq. (|4Tj) we obtain 



FIG. 7. Plot of [f(e - hw) - f(e + hw)]/(huj) (solid line) 
and —df(e)/de (dashed line) as a function of e — 6f for 
T = 0.43 meV, hw = 1.17 meV. 




FIG. 8. Plot of the convolution / p(e - Kb/2) p(e + MJ/2) de 
(solid line) and of the Lorentzian profile (P /-kK)/(£j 2 + F 2 ) 
(dashed line) as a function of Kj/"/b for P = (37r/8)(7s/ft). 



Re aiJ (cj) = 2A/Ve 2 / de 



f(e-tiu)/2)- f(e + Hu/2) 
hw 

x h p(e — hw/2 — e{) p(e + haj/2 — e;/) 
1,1' 



x^2\(l,a\v x \l',a')\' 

a,ot' 



(56) 



In this expression the disorder enters only in the den- 
sities of states, p(e), while the velosity matrix elements 
are the same as in the absence of disorder. Moreover, 
when the quasiclassical approximation is valid, (e; de- 
pends smoothly on I) they can be calculated classically. 
Indeed, comparing Eqs. (pTtj) and ((55)) . one can establish 
the correspondence 



df ;; f(e-Hu/2)-f(e + hw/2) 
de hcu 

Ll X Lly 



T/ir 

) 2 + r 2 



<-> hi p(e) p(e + hw) de, 



(57a) 
(57b) 

(57c) 



where we denoted Cj — uj — nuj c or w — (e;+„ — e{)/%. 
The hrst two lines hold reasonably well in our case (even 
though T = 5 K = 0.43 meV < hw = 1.17 meV, see 
Fig. [7]), while in the third one the convolution of two 
semicircles is somewhat different from a Lorentzian, so 
there is no unique relation between Y and 7b- We choose 
to match the peak heights (that is, the value at Cj = 0), 
which fixes 



37T j B 

II' 



(58) 



as illustrated by Fig. [5] 



only the first term of the ladder series in Fig. [^c) con- 
tributes to the inter-Landau-level transition, in full anal- 
ogy with the case of the parabolic spectrum 0]. 
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